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1 . Introduction 


1.1 Overview 

The prediction of the behavior of turbulent reacting flows 
la a major concern of current combuption research. Such modeling 
is required in order to enhance the understanding of the phenomena 
involved and to design and evaluate the performance of combustion 
devices, A principal element to be derived from the modeling 
of these flows is an expression for the reaction rates of the 
various species ln.volved in any particular combustion process 
under consideration. 

Currently, several approaches for the determination of 
the properties of turbulent, reacting flows exist. Of the 
present approaches, the method of utilizing an assumed probability 
density function (pdf) for temperature and species, concentrations 
has been selected for use in this study. The motivation for 
selecting this approach is its relative computational simplicity 
and its basis in the probabilistic nature of turbulence, 

1.2 Objective 

* 

This work is examining the effects of temperature and 
species concentrations fluctuations on reaction rates in turbulent, 
reacting flows by means of the presumed pdf approach. The follow- 
ing items are the subject of the present study: 

1. The utilization of the temperature-derived most-likely 
pdf of reference (2) to describe the effects of nampera- 
ture fluctuations on the Arrhenius reaction rate 


constant. 


, 2. The utilization of the most-likely bivariate pdf 

of reference (2) to describe effects of temperature 

» 

and species concentrations fluctuations on the 
reaction rate. 

3. Development of a criterion for the use of an 
"appropriate" temperature pdf. 

Section 2 presents a brief review of the general theory 
of pdf's. Sections 3 and 4 review the formulation of models 
to calculate the mean turbulent Arrhenius reaction rate constant 
and the mean turbulent reaction rate, respectively. Section 5 
presents the results of calculations using these models. 

Section 6 presents conclusions and the direction of future work, 
dealing with the combined effects of temperature and species 
fluctuations on reaction rates. In Section 7, we deal with the 
criterion for the selection of a temperature pdf. 


2. Review of the General Theory 
The discussion in this section is taken largely from 
reference (1) . 


2.1 Probability Density Functions 

A pareuneter x is said to be a continuous random variable if 
there exists a probability density function, p(x), which satisfied 


the following conditions. 


(31 


JsCx) i O -for ad X 


(la) 
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The pdf may, of course, be defined on an interval other 
than (-«, + «) , for example (0,1). Since the pdf is defined on 
a specific interval, the functional value of p(x) is zero 
elsewhere. 

Equation (1) must be satisfied by a pdf of a one-dimensional 
continuous random variable. Probability density functions may 
be written for a multi-dimensional continuous random variable. 

A probability density function for a two-dimensional continuous 
random variable, denoted pCx,y) , 'is termed a bivariate or joint 
pdf. For such a pdf, the conditions corrvesponding to equation 
(1) are: 


> o 

for all X and y 

/ f * 1 


(2a) 


(2b) 


~|0 

The expected, or mean value of a one-dimensional continuous 
random variable, x, is expressed as: 




(3) 


The variance of a one-dimensional continuous random variable, 
X, is expressed as: 

■o 

= VCO = r 

( 4 ) 
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. A useful quantity derived from the variance ia the standard 

* i!:: . 

deviation which is the square root of the variance. Equation 
(4) is a special case of the more general expression for th^ 

Kth moment of a random variable, x, about its expectation. 

The general expression for the Kth moment of a random variable 
is given by : 

* f( X -EU)) • l>U) 

Clearly, for the special case of K=2, equation (5) is equal 
to equation (4) . 

Equation C3) may be extended to functions of a continuous 
random variable. In the case of a function f(x) of a one- 
dimensional random variable having a pdf P(x) , the mean value of 
f(x) is: ^ 

fcO * f * f> ^ ^ 

similarly, in the case of a function g(x,y) of a tv/o-dimensional 

* 

random variable, having the joint pdf ptx,y), the mean value of 
g(x,y) is: - ^ 

.«o 

The correlation coefficient, is a parameter defined 

with respect to the two-dimensional/continuous random variable 
(x,y) and is expressed as: 


( 8 ) 


P 

• iky 


/V<x) N/ty) ’ 


The correlation coefficient ia a measure of the degree of 
linearity between x and y. Values of the correlation co- i 

efficient near +1 or -1 reflect a high degree of linearity, 
while values of the correlation coefficient near ssero indicate 
a lack of linearity. Positive values of the correlation co- 
efficient indicate that as y increases, x increases. Negative 
values of the correlation coefficient indicate that y increases 
as X decreases. 

The numerator of equation (8) is defined as the covariance 
of X and y. The covariance is denoted by and is expressed as: 






, . C9) 

The significance of the covariance can be ascertained by con- 
sidering a two-dimensional random variable (x,y) . x and y are 
termed independent random variables if the value of x has no 
influence on the value of y (and likewise, the value of y has no 
influence on the value of x) . When x and y are independent 
random variables the covariance is zero. Hence, the covariance 
may be considered as a minimum ’'criterion" of statistical 
dependence. This criterion can assure, at the very least, that 
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vhen the covarience is .not zero, the variables are not indepen- 
dent. However, no statement can be made concerning independence 
on this basis alone, i£ the covariance is zero. 


2.2 Arrhenius Reaction Rate Constant and Reaction Rate 
Expressions 

For the purposes of the study presented herein, the follow- 
ing one-step, irreversible reaction is considered: 

r + ^ o -•‘p 

( 10 ) 


F, 0, P *= fuel oxidizer, and product species 
respectively 

0 - stoichiometric coefficient 
The disappearance of P and O are related by: 




JL iii 
O At 


( 11 ) 


where: c„, °n “ concentrations of fuel and oxidizer, 

respectively 

The negative signs in equation (11) indicate that the concentra 
tions of P and 0 decrease as the reaction proceeds. The 
reaction rate of fuel is expressed as: 



At 



* t 

where k(T) = Arrhenius reaction rate constant (12) 

T = temperature 

n, m = constants dependent upon the particular reaction. 

The Arrhenius reaction rate constant is expressed as: 

k(T) = AT®exp (-T^/T) (13) 


6 - 


where A, B « constants dependent upon the particular reaction 
* activation temperature of the particular reaction. 
Equation (13) can be modified by the introduction of the dimen- 
sionless temperature: 


■fc »s 


T - Tmin 


Tmax - Tmin 

M , 

The temperatures Tmax and Tmin are defined so that the 
dimensionless temperature t can only assume values within the 
interval (0,1). Thus, Tmin may be the lowest temperature of 
the unreacted constituents and Tmax may be the equilibrium 
combustion temperature of the reaction. Substitution of 
(14) into (13) yields: 

KCt) •* A Tk /(K, t 4 


(14) 


where 


Tmax - Tmin 


^2 = Tmin 


Equation (12) is also modified by the introduction of the 
following dimensionless concentrations: 


(15) 




Cp. -Cp. 


■=* 

I o 


c - c 

’“o “ 


In these expressions, ''max'' and "min" denote maximum and 


(16) 


minimum values, respectively. Equation (16) ensures that and 
r“Q only assume values within the interval (0,1). For the 
case where the minimum concentrations are zero and the maximum 
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concentrations are the initial values, the combination of 
equations (12) and (16) yields: 

Cr.c:^ 

(17) 

* . 

where c„® Cq® « initial concentrations of fuel and oxidizer, 

' respectively 

2.3 Mean Turbulent Arrhenius Reaction Rate Constant and 
Mean Turbulent Reaction Rate Expressions 

In a turbulent reacting flow, the mean turbulent Arrhenius 

reaction rate constant can be calculated by treating the 

temperature c,s a continuous random variable and specifying an 

appropriate probability density function for the temperature. 

This concept is employed in Reference (1) along with the 

analogs of equations (6) and (15) in the present work to yield: 

f 

K.U) . r K ttv fCO Jt 

o 

(18) 

where ; = mean turbulent Arrhenius reaction 

rate constant. , ... 

The above expression provides a direct, relatively simple 

method of taking into account the effects of temperature 

fluctuations on the Arrhenius reaction rate constant in a 

turbulent, reacting flow. 

Similarly, the expression for the mean turbulent reaction 
rate can be developed from an extension and combination of 
equations (7) and (17). Thus: 

• I • 

/ f f '^F p 

® o O ygj 
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where Wp ■ mean turbulent reaction rate of fuel. 

This expression utilizes a joint pdf for temperature and species. 
Equations (18) and (19) can be applied to one-step reactions, 
or multi-step mechanisms, by considering each elementary reaction 
separately. The pdf's used here are considered to be valid at 
an instant in time, and thus are not functions of time. 

Xn order to compare the magnitude of a mean turbulent 
Arrhenius reaction rate constant to the corresponding laminar 
term, an amplification ratio is defined. The laminar Arrhenius 
reaction rate constant is calculated by inserting the mean 
dimensionless temperature in equation (15). The amplification 
ratio is obtained by dividing equation (18) by (15) ; 


Z « 

where? z m reaetien rate constant amplification ratio’ 

Jc^ = mean turbulent reaction rate constant 

«= corresponding "laminar" reaction rate constant 


( 20 ) 


A similar term may be defined for the ratio of a mean turbulent 
reaction rate to the corresponding "laminar" reaction rate. The 
"laminar" reaction rate is calculat<ad by inserting mean values 
of dimensionless temperature and concentrations into equation., . (17) 

The ratio of equation (19) to (17) defines the reaction rate 
amplification ratio: 

z\ - w / (21) 

I 

vhere Z = reaction rate amplification ratio 
w^ *= mean turbulent reaction rate 


Wjj, = corresponding "laminar" reaction rate 


Equation (21) is the reaction rate amplification ratio 
in which the mean turbulent reaction rate is calcuilated with 
consideration of the combined effects of temperature and species 
concentrations? fluctuations, in order to compare this amplifi- 
cation ratio to one which considers only temperature fluctuations , 
the following is defined with the aid of equation (17) : 

) Cr, Co) 

This expression involves a mean turbulent reaction rate in which 
a mean turbulent reacton rate constant is employed along with 
laminar values of concentrations. The corresponding laminar 
reaction rate also contains these same concentration terms.. 

Thus, the concentration terms cancel. The reaction rate 
amplification ratio for this case is: 



( 22 ) 

This amplification ratio may be compared to that calculated . . 
by use of equation (21). Hence, the combined effects of tempera- 
ture and species concentrations fluctuations on the reaction 
rate may be compared to effects of only temperature fluctuations. 

3. Probability Density Functions 
for Temperature 

3.1 Most - Likely pdf 

The most-likely pdf is shown in reference [2] to be the 
statistically-most-likely pdf of a continuous random variable 
for a given flow. It is important to note that although this 
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pdf is the most-likely pdf, it is not necessarily the pdf 
for a given flow. The utility of the most-likely pdf is 
that any number of moments can be incorporated into the pdf 
and, hence* greater accuracy achieved. 

Two possible shaxjes of the most-likely pdf for one 
variable (e.g., temperature) are shown in Figure 1 

P<t) 


OL 
0 

i 

Figure 1. Possible shapes of the one-variable 

most-likely pdf (a) spike (b) Gaussian-like. 



1.0 



The expression for the temperature-derived most-likely pdf is: 

Iptt^ s j X,t ] 

* 

C23) 

where the constant coefficients lO, 11, 1?. are obtained 

from the simultaneous solution of the following constraint 
equations for known values of dimens ion- less mean temperature 
and the mean square fluctuation of dimensionless temperature; 
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I 

[ t» d-t • \ 

o 

[ fc « t 

o 

* X 

f U - 1 ^ • (>CO s t 


( 24 ) 


(25) 


(26) 

Equation (23) may be utilized in equation (18) to obtain 
the following expression for the mean turbulent reaction rate 
constant: ^ 

Kt ’ A [ (K.t ■ »»(> [-T». I (K,t ♦Kv'i]- |> <l't 

* C271 

where pCt) is the temperature-derived most-lijcely pdf, 

, ••.equation (23). 

The corresponding "laminar" term is calculated by use of t 
in equation (15). This yields: 

^ 4-K I 

The ratio of (27) to (28) is the reaction rate amplification 
ratio, Z: 

(29) 

For the purposes of a parameter study, values of mean 
temperature and mean square temperature fluctuation are selected 
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* in ord@K to ©jcamino th© effocts oU thoir variation on the mean 
turbulent Arrhenius reaction rate constant, 'J?he results of 
this study are presented in Section 5. 

The computational procedure to determine the values of 
the reaction rate constant amplification ratio is described 
below; 

Procedure ; 

~ " 2 

1. Select values for A, B, TA, Train, Tmax^ t and t' 

2. utilize Newton's Method to solve for the constant 
coefficients of equation (23) . 

3. Numerically integrate equation (27) to obtain 

4. Utilize t in equation (20) to obtain IcA 

5. Calculate Z with equation (29). 

3.2 Beta pdf 

The beta pdf is utilized in Reference (1) as a pdf for the 
temperature. The selection of the beta pdf is based on its 
usefulness in approximating experimentally determined temperature 
pdf's in turbulent? non-reacting flows. Two possible sheipes 
bf the beta pdf are shown in Figure 2, 

1 . ‘ ■■ 

I 



Tlie expression for the beta pdf is: 


l>CO 


rcaVPCfe^ 


t^^Ci-ty ; 04 til 


-fc-i 


vhere P » gairana function 

— * *2 ■ 

a, b « constants dependent on t and t' 
The expressions for OL and b are: 

a- t[tO-tVt^->l 

b. 0-o[t0^t)/v‘-i] 


( 30 ) 


Ola) 


(31b) 

Equation (30 is utilized in equation (18) to obtain 
the following expression for the mean turbulent Arrhenius 

reaction rate constant: 

I 

* ■* t32) 

where p(t) is the beta pdf. 

The procedure for evaluating equation (32) is given in 
Reference (1) . Equation (32) and equation (28) are inserted 
in Equation (29) to obtain the reaction rate constant amplifi-* 
cation ratio. 

The results of a parameter study with the model given by 
equation (32) are presented in Reference (1) . In this para- 
meter study, values of mean temperature and mean squaro 
temperature fluctuation are selected in order to examine effects 
of their variation on the mean turbulent Arrhenius reaction 
rate constant. The results of this study are presented in 


Section 5 for comparison with the resuXts Obtained by utilizing 
the temperature-derived most-likeXy pdf in equation (29) . 


4. iJoint Probability Density Functions for 
Temperature and Species 


4.1 Most-likely Bivariate pdf 

The most-likely pdf of Section 3 is easily extended to 
two or more variables (e.g., temperature and one species 
concentration; temperature and two species concentrations, etc.). 
The form of the most-likely bivariate pdf for temperature and 


one species concentration is: 


.( 2 ) 




Where <1 «» a priori probability dependent upon the' 

reaction rate 

t *■ dimensionless temperature 

« dimensionless concentration of fuel, II 2 
constants 

The constant coefficients V® ,Xi determined 

from the simultaneous solution of the foll owing constraint 


equations for known values of ^ 

I I 

) [ t> C«-p .t) jTp at • 1 


0 o 

1 f 


i ^ 'V • |»<'f ‘‘•V 


O o 
• I 

o o 




(33) 


(34) 


(35) 


( 36 ) 


•IS- 


.( 37 ) 






Rguation (34) follows from equation t2) , Equations (35) and 
(36) are expressions for the dimensionless mean concentration 
of fuel species and dimensionless mean temperature, respectively. 
Equation (37) is the expression for the covariance of the 
dimensionless fuel species concentration Tp and the dimensionless 
temperature t . 

For a given value of t and fp , the value of the covariance 
is obtained by utilizing equation (8). Feir the purposes of 
this study, an approximate value of the correlation is taken 
to be -0.9 from Reference (1), for an assumed adiabatic, 
turbulent diffusion flame . Hence , 






V 






(38) 


where 






S3 -0.9 


In addition, in order to utilize equations (35) , (36) , (37) ,' 
and (38), the following relationships from Reference (1) arc 
employed : 

ri» « 


r-F • T ^ 


(39) 


( 40 ) 


3?he term appearing in equation (33) is an a priori 

probability dependent upon the reaction rate. In the present 


study, t and rp are treated as passive scalars. This 

/)' 

simplification is seen as an initial step in the utilization 
of the most-likely pdf. For this case, Reference (2) shows 
that the q-term is a constant. For purposes of calculation, 
this term is set equal to unity. 

Equation C33) may be utilized in equation (19) to obtain 
the following expression for the me/an turbulent reaction rate: 

— * * ^ r *1 

^ H VTK/(K,t*ic^)Jr^<:,* . 

e ft 

[jc; crp-o+ c* V 

( 41 ) 

where p(t,rp) = most-likely bivariate pdf. 

The corresponding "laminar” term is calculated by use of t and 
fp in equation (17). Refeeense (1) shows that equation (17) 

may be rewritten as: 

( 42 ) 

The ratio of equation (41) to (42) is the reaction rate 

I 

amplification ratio Z : 

7 * . 

hm ^ j 

w ^ (21) 

For the pu)irposes of a parameter study with the model 
given by equation (41) , the following one-step, irreversible 
reaction is considered: 

(43) 
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The jresults of this parameter study are presented in Section 5. 

The computational procedure to determine the values of 
the reaction rate amplification ratio is described below. 

I 

Procedure 

1. Choose values for A, B, T^, Tmini Tmax# t and t'^ 

2. Utilize Hewton's Method to solve for the constant 

coefficients of equation (33) , ^ 

3. Numerically integrate equation (41) to obtain w^. 

4. Insert t and r into equation (42) to obtain w, . 

5. Calculate Z with equation (21). 


4.2 Modified Bivariate Gaussian pdf 

This is an alternative two-dimensional pdf, as initially 
presented in Reference (1) and is given by equation (44) : 


■f C-t.r'i » -r-; 


o 4 -fe s 

ot I 


(44) 


In the above equation, p(t,r) is the bivariate Gaussian pdf 
given by: 


I? Ct.r') = ; ; 

an «;<r^ 



■ - "V 





- OO <, t < pO 
-.eo 4, *- 4 oo 
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Eqiiatlon (44) is utilized in equation (19) to obtain 




the following expression for the mean turbulent reaction rate; 

I I 


0 6 ^ 


( 45 ) 


where (t,rp) « modified bivariate Gaussian pdf. 

Discussion of the modified bivfrriate Gaussian pdf and 
the procedure for utilizing equation (45) are given in 
Reference (1) . 

The ratio of equation (45) to (42) is the reaction rate 

t 

amplification ratio, Z : 

z*“ 


w 


(21) 


The results of a parameter study with this model are 
presented in Reference (1) . This parameter study considers 
the one-step, irreversible reaction given by equation (43) . 

The results o| the parameter study of Reference (1) are 
compared to those obtained by means of the most-likely bivariate 
pdf in Section 5. 


5, Results 


Throughout the discussion in this section, the use of 
the terms "temperature” and "concentration" denote mean dimen- 
sionless temperature and mean dimensionless species concentra- 
tion, respectively. The use of "r" in all figusies denotes 
r. 


•H, 




- 19 - 


. u 




iripieiafiiri» 


5il The Effects of Temperature Fluctuations on the Arrhenius 
Reaction Rate Constant 

The effects of temperature fluctuations on the Arrhenius 
reaction rate constant are presented in this first subsection. 
The calculations are performed with the temperature-derived 
most-likely pdf. The results of these calculations are com- 
pared with those obtained by use of the beta pdf taken 
from Reference (1) . 

5.1.1 Computational Aspects 

All numerical integrations are performed using Simpson's 
Rule. Newton's Method is used to solve for the unknown 
co'^efflcients in equation (23) . Details concerning this 
technique and the associated proof of convergence and error 
analyses will appear in the next progress report. 

5.1.2 Discussion of Results 

Table 1. shows the cases considered. The data presented 
are obtained from Reference (5) and are considered to be typical. 
All the following results are stated with regard to both the 
temperature-derived most-likely and beta pdf's. Proposed ex- 
planations for the phenomena involved are given in Ref<;rence (1) . 


Table 1. Case Data 


Case 

No. 

B 

^A 

■'max 

(°K) 

1 

0 

10116. 

2500. 

2 

0 

5000. 

2500. 

3 

1 / 

10116 . 

2500. 

4 

& 

10116. 

3000. 


Note: A = 8.4 X 10^^, T . = 500°K, for all cases. 

* TTl T * 


Vairiatlon with mean temperature and temperature fluctuation 

Figures 3. through 7. show the variation of the reaction 
rate constant amplification ratio with mean temperature and 
mean square temperature fluctuation. The following are 
observed : 

1) The reaction rate constant amplification ratio 
increases with increasing values of mean square 
temperature fluctuation, at constant values of mean 
temperature for both pdf's. Refer to Figures 3. 
through 7. 

2) The reaction rate constant amplification ratio 
decreases with increasing values of mean tempera- 
ture fluctuation for both pdf's. Refer to 
Figures 3. through 7. 

3) The same values of Z are obtained with the use of 
both pdf's for values of t « 0.5, 0.6. Refer to 
Figures 5. and 6. 

4) Nearly the same values of Z are obtained v;ith the 

use of both pdf's for values of t equal to 0.2 

- ~~2 

through 0.7 and, for t = 0.0 and t* £ 0.04. 

Variation with activation temperature 

Figure 8. shows the variation of the reaction rate constant 
amplification ratio with the activation temperature. As may be 
seen, a decrease in the activation temperature results in a 
similar decrease in Z for the tvro pdf's. 
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Variation with tv^mperatiire exponents B . 

Figure 9. shows the variation of the reaction rate 
constant zunplification ratio with the temperature exponent, 

It is observed that an increase in the temperature exponent 
results in similar increases in A, for both pdf's. 

Variation with maximum temperature 

Figures 10, and 11. show the variation of the reaction 
rate constant amplification ratio with the maximum temperature. 
It is seen that in both cases the value of Z is relatively 
insensitive to changes in Tmax. 

Prom these results, it is concluded that quite similar 
results are obtained by the use of the beta and most-likely pdf' 
Since the former has been experimentally verified in many non- 
reacting, turbulent flows, this similarity in behavior lends 
support to the physical correctness of the most-likely pdf, 

5.2 Effects of 'temperature and Species Fluctuations on the 

Reaction Rate 

The effects of temperature and species concentrations 
fluctuations on the mean turbulent reaction rate are presented 
in this subsection. The calculations are performed with the 
most-likely bivariate pdf. ' ' The results of these calculations 
are compared with those obtained by means of the modified 
bivariate Gaussian pdf take;}! from Reference (1) . 

5.2.1 Computational Aspects 

During the course of this study, a method was developed to 
transform the double integrals of equations (34) through (37) 
to single integrals. This rendered the evaluation of the asso- 
ciaiied unkno^m constants practical by means similar to those 
mentioned in Section 5.1.1. Details of this transform method 





and the associated numerical techniques will appear in a future 
status report* 

5.2,2 Discussion of Results 

Table 2. shows the cases considered. The initial values 
of concentrations are taken from Reference (5) , All other data 
are also taken from Reference (5) and are considered to be 
typical. All the following results are stated with regard to 
both the most-likely bivariate and modified bivariate Gaussian 
pdf's. Proposed explanations for the pheonomena involved are 
given in Reference (1) . 


Case 
, .No. 

B 

^A 

(°K) 

^max 

(°K) 

P 

»2 

«2 

t ' 
vs 

2 

1 

0 

f. f 

‘10116. 

2500. 

-0.9 

«2 

« 2/3 

f' 

4 

0 

10116 . 

3000. 

-0.9 


It 


5 

0 

10116. 

2500. 

o 

• 

CO 


M 


6 

0 

10116 . 

2500. 

-0.9 

r'^ 

«2 

« 1/2 

t' 


Note: A « 8.4 x 10^^* Train = 500°K, C° = 2.405 x IQ-^gmol/cm^ , 

“2 

Cq = 1.5 X 10”®gmol/cm^ , = 1-t, for all cases. 


Variation with fluctuations and mean values of temperature and species 



! 


I 



i 


Figures 12. through 17. show the variation of the reaction 
rate amplification ratio with temperature and species fluctuations. 


at .various constant values of mean temperature and mean concentra- 
tion. l?hese results consider combined effects of temperature and 
species fluctuations, tPhe following are observed; 

1) The reaction rate amplification ratio increases vfith 
Increasing values of fluctuations Jfor values of t from 
C.2 through 0.5. Refer to Figures 12. through 15, 

2) The reaction rate amplification ratio reaches a 
relative maximum and then decreases with increasing 
fluctuations at t «• 0,6. Refer to Figure 16. 

3) The reaction rate amplification ratio decreases with 
increasing fluctuations and is l ^/ss than unity at 

t » 0.7. Refer to Figure 17 . implication of Z < ^ 

is that species fluctuations are resulting in an 
"unmixedness" effect which is reducing the turbulent 
reaction rate to values less than the “laminar ''-calculated 
values . 

4 ) The reaction rate amplification ratio decreases with 
Increasing values of t Cand correspondingly decreasing 
values of r) at constant values of fluctuations. Refer 
to Figures 12. through 17. This implies that in this 
turbulent flame, the principal turbulent amplification 
is occurring in the region of flame initiation. 

5) The use of the most-likely bivariate pdf yields values 

I 

of Z greater than, equal to and less than values obtained 
with the use of the modified bivariate Gaussian pdf, for 
values of t “ 0.2, 0.3 and 0.4, Refer to Figures 
12. , 13. , and 14 . 
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. 6i The use of the most-likely bivariate pdt always yields 

I 

values of Z greater than the values obtained with 
the use of the modified bivariate Gaussian pdf for 
values of t 0.5, 0.6, 0.7. Refer to Figures 15., 

16. and 17. The results obtained in the previous 
section lend greater credence to thta values obtained 
using the most-likely bivariate pdf. 

Figure 18 reveals that both joint pdf's are sensitive to 
the value selected fo,r which v’as not the case for the 

temperature-only pdf ' s . 

Variation with the correlation coefficient 

Figure 19, shows the variation of the reaction rate 
amplification ratio with the correlation coefficient. It is 
observed that the value of Z* exhibits a relatively small decrease 
with an 11% increase in the >7alue of y£>, 

~~ 72 * , 2 

Variation with r „ vs t 
^ H2 ■ 

Figure 20. shows the variation in the reaction rate ampli- 
fication ratio with an approximate change of 30% (2/3 to 1/2) , in 

. .r _,2 

the relationship between r h 2 and t . The following are observed 

I 

1) The value of Z is generally lowered by this. 

I 

2) The decrease in Z is substantial with the use of the 
most-likely bivariate pdf. This suggests that the most- 
likely pdf will serve as a more sensitive indicator of 
the correctness of a given turbulence-reaction model, 
when compared with the modified Gaussian. 
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5 ..3.1 Discussion of Results 


Variation with fluctuations and mean temperature 

Figures 21. through 26. show coraparisons of the variation 
of the reaction rate amplification ratio for the two-variable 
and. one-variable models with values of the fluctuations and 
mean temperature . 

The following are observed: 

1) The reaction rate amplification ratio increases with 
increasing temperature fluctuations at constant values 
of t, for t = 0.2 through 0.5. This trend is realized 
for both the one-variable and two-variable models. 

Refer to Figures 21. through 24. 

2) The value of Z decreases with increasing t at constant 

■ — - A - 

values of t' for both the one-variable and two-variable 
models. Refer to Figures 21, through 26. 

I 

3) The value of Z increases with increasing temperature 
fluctuations at t = 0.6, with the use of the one-variable 

^ f 

model. Also at this value of t, Z increases, reaches 
a relative maximum, then decreases with increasing 
temperature fluctuations, with the use of the two-variable 
model. Refer to Figure 25. 

I 

4) The value of Z increases with increasing temperature 
fluctuations at t = 0.7, with the use of the one-variable 

' model. For all values of t'^, this model yields values 
of Z > 1. At this same value of t, z' decreases with 
increasing temperature fluctuations and is alv/ays less 
than unity, with the use of the two-variable model. 

Refer to Figure 26. 
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These results clearly indicate the necessity for multi- 
dimensional pdf models since, in temperature-only models, 
the species "unmixedness” reduction in the reaction rates 
(2 < 1) is totally absent. 

6. Conclusions and D irect ion of Future Work 

• mm, I l l- .III 

6 . 1 Conclusions 

6.1.1 Effect of Temperature Fluctuations on the Arrhenius 
Reaction Rate Constant 

The effects of temperature fluctuations on the Arrhenius 

I 

reaction rate constant are assessed by treating the temperature 
as a continuous random variable and utilizing two pdf's for 
temperature. The pdf's are the temperature-derived most-likely 
and the beta. The ratio of the mean turbulent Arrhenius res,ction 
rate constant to the corresponding laminar term is defined as the 
amplification ratio, 2. The results obtained by means of both 
pdf's exhibit the following trends: 

1. 2 is always greater than unity.’ 

2. 2 increases as t'^ increases, at constant t 

3. 2 decreases as.t increases, at constant t"’"^' 

4. For wide ranges of t and t , the use of both pdf's 
yields nearly identical numerical results. 

5. 2 is relatively insensitive to moderate changes in T« 
and B. 

6. 2 is insensitive to moderate changes in 

6.1.2 Effects of Temperature and Species Fluctuations on the 
Reaction Rate 

The effects of temperature and species fluctuations on 
the reaction rate are assessed by treating both the temperature 


and species concentration as continuous random variables 
and utilizing joint pdf’s relating them. The pdf’fs are the 
most-likely bivariate and modified bivariate Gaus/sian. 

A reaction rate amplification ratio zVis defined at the 
ratio of a mean turbulent reaction rate to the corresponding 
laminar term. Also defined is a reaction rate amplification 
ratio with a mean turbulent reaction rate composed of a mean 
turbulent Arrhenius reaction rate constant and laminar values 
of concentration. This amplification ratio only considers the 
effect of temperature fluctuations on the reaction rate. It 
is possible to compare the combined effects of temperature 

f 

and species fluctuations to those of only temperature fluctua- 
tions by comparison of the two-variable (temperature and species) 
and the one-variable (temperature) models. These one-variable 
models utilize the temperature derived most-likely and beta 
pdf's. The results of the two-variable models exhibit the 
following trends: 

I 

1. Z of the combined case may be greater than, equal to 

or less than that of the temperature-only case. 

• * ^ 

2. Z of the combined case increases with increasing 

temperature and species concentration fluctuations, at 

low values of mean ten^perature and correspondingly high 

values of mean species concentration. At high values 

of temperature and correspondingly low values of con- 
» 

centration, Z decreases with increasing fluctuations. 

I 

3. Z of the combined case is sensitive to changes in T 

^ max. 

I 

4. Z of the combined case is nearly insensitive to small 
changes in p. 
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I 

h, The values of z obtained with the use of the wost- 
likely bivariate pdf are moderately sensitive to 
changes in the relationship between temperature and 
species fluctuations. The values of Z* obtained with 
the use of the modified bivariat^ Gaussian pdf is in~ 
sensitive to moderate changes in this same relationship, 

6,2 Direction of Future Work 

1. Continue development of a temperature-two species 
fluctuations model (previously identified as model II, 

t 

Reference Ic) . 

2. Evaluate the "q" term in equation (33) to allow the 
treatment of reactive scalars. 

3. Utilize the three-variable analog of equation (33) in 
order to eliminate the assumption on which equatibn (41) 
is based (previously identified as model HI, 

Reference Ic) . 

I 

7. Development of a Criterion for the use of 
the 'l^'amp''‘R'bf ' ' 

♦ 6 . 

7 . 1 Introduction 

A criterion for selecting betv>reen the beta (or oxie-dimensiq/nal 
most-likely) and ramp pdf's is under investigation. The flatness 
factor and skewness are studied for potential use as criteria. 

The ramp pdf generates values of flatness factor between 2.0 and 
2.5. The beta pdf also generates these values of the flatness 
factor since it is a continuous distribution (its range of flat- 
ness is 1.66 and 8.64.) However » when the flatness factor is 
between 2.0 and 2.5, the skewness for each of the pdf's is 


aistinot. one to chooso sHownoss as a sopondaKy 

oKitQKion vith flatnoss factor boin^r prib^ary. 

A coraparisoa of the beta and wost-liheXy pdf's is made with 
the conclusion that they do not yield equivalent valvies for the 
higher order moments. 

7.2 Moments of a pdf 

Every probability density function (pdf) has an infinite 
number of moments associated with it. The a\ost coiwt\only used 
moment is the first, called the mean, vi. The mean is defined 


yiL. f 


Alternatively, u •» t may be used. This moment is also 
called the expected value, E(t). 

The second moment is the variance, V. The square root of 

the variance is the standard deviation, u, and is given by: 

#• 

(t‘-» 

If the polynomial is expanded and multiplied through b5>* pCt)', 
one obtains 

(T'* • ^ f>Ct) 4'^ 

The only factor that is not known in Equation (18) is 

•« 

£ (tO - ( 



similarly, higher order moments will also contain a term of 
the same form» or 


£ (t*) X 


( 50 ) 


where V: is> the moment desired. 

In order to calculate any moment, it is sufficient to 
calculate B(t ) and then multiply by the appropriate constants 
to obtain the moment desired. 

Moments higher than 2 are given by the following general 
equation : 

•* ■ M 

Note that expansion of the (t-n)**^ term will yield the term 
t** p^^^ dt. 

*fhe third memenHi is called the skewness, and the fourth 
is called the flatness factor. 


(51) 


7.3 Moments of the Beta pdf 

The beta pdf, proposed by Rhodes (6) to describe turbulent 
flow properties, is given by: 

V 


where 


and 


^(t)xj:i o-t) 


dt 


O 


(52) 


a» t 


f ^ 0 - i ) 7 

L^-'J 


(53) 


L ?= 7 
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The denominator of Equation (52) can be expressed using 
the Gamma function; 

! i., . I'M- ret) 

i ‘ <'-'1 = — i .. o 


( 55 ) 


The moments of the Beta pdf are 

« 


o 

where V is the mean. If the polynomial in Equation (56) is 
expanded, there is always a / t^ p^^^ dt term to be evaluated. 
It can be evaluated using a form of Equation (55) ; 


(56) 


ECi*) 


r^+fc) 


/ 1 J-t 


P(«) -PCbi 


n Cet-b) 




r<*)- p(fc> 




(57) 


tS8) 


If one lets x in Equation (57) equal (k+a) and y«b, then 

/ . it\ ^ p • n(c) 

^ ^ *P(a).PC(») ’ 

^ r P(e). P (k ) 


( 59 ) 
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Equation C59) has been used to evaluate the moments. 

The skewness and flatness factor distributions as functions 

■ y 

of t and t' are shown in Figures 27 - 30. The distribution 
is symmetric about t»0.5, therefore only one graph is needed 
for both t=0.4 and t«0.6 since they have the same values. The 
only difference is in the skewness ~ for t=0.4, the skewness 
is positive, for t=0.6, the skewness has the same magnitude but 
is opposite in sign. The same holds true for t«0.1,0.9; t=0.2,0.8; 
t«0.3,'0.'7, etc. 


7.4 Moments of the Pope Pdf 

The most-likely pdf proposed by Pope^^'^^ h-fis the form; 

pCi)* t 

where ^1^ ”2 three constants whxch are functxons of 
t and t' . Since there are three unknown^, three equations are 
required to determine the three constants. However, only 
two algebraic equations can be developed (Appendix B) : 

tM'l 

_ I 


(60) 


(61) 


(62) 


The required third equation is an integral equation defining 
one of the properties of a pdf, or its first two moments, namely 

t 




( 63 ) 


( 64 ) 


|t- i 


I 

Jct-i*)” .f,ai 

o ' 

(65) 

Any one of these three equations can be used as the third 
equation. 


Once the three X's are known, the higher order moments 


can be found using the following equation: 

I 




0 ( 66 ) 

Graphs of the skevmess and flatness factor as a function 

» ■■ W 

of t and t' are presented in Figures 27 - 30, As with the beta 
pdf, the values of the flatness factor are symmetric about 
t=0.5. That is, the same values of the flatness factor are 
obtained for t=0.4 and t=0.6. The same is true for the skew- 


ness, except that the signs are reversed (negative values when 


'b=0.5). One interesting phenomenon encountered was that there 
seemed to be no solution obtainable for some t when t=0.1 or 
0.9. Some data points were obtained, but not enough to permit 
a comparison with other pdf's. 


7.5 Comparison between the Beta vs. Pope Pdf's 

The ,l!<^yior of the beta and Pope pdf's is very similar 

near t=0.5. As can be seen from Figure 30, both curves follow 

2 - 

the same trend and approach each other at high t' . At t=0.4,0.6. 
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the flatness factor curves are fairly similar and approach 
each other at high values of the fluctuations. The skewness 
curves are not as similar as the flatness curves. 

The differences between the pdf's becomes more pronounced 
at t«0.3,Q.7. The skewness curves are essentially the same 
except they are a little further apart. The flatness factor 
curves are distinct now as the Pope pdf has a maximum while the 
beta does not. At the maximum, the flatness factor for the 
Pope pdf is approximately 13% higher than for the beta pdf. 

At t=0.2,0.8, t;here is a large difference in the flatness 
factor curves. Both curves have a maximum, but the Pope pdf 
has a much more pronounced peak. The skewness curves have 
basically the same shape as before, but again are further apart. 

The two skevmess curves always cross each other, with the 
point of intersection moving to the right with increasing t. 
After the intersection, the Pope pdf yields higher values than 
the beta pdf. 

The conclusion that can be drawn from these curves is that 
the beta and Pope most-likely pdf's yield values for the higher- 
order moments which can differ substantially. Both encompass 
the same range of values, but the differences are such that, on 
this basis, they are not interchangeable, 

7.6 Moments of the Ramp pdf 

The ramp pdf was developed by Antonia and Atkinson after 
observing ramp-like temperature fluctuations under a variety of 
conditions. The pdf is composed of two parts - a raraplike 
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structure with high -freguency Gaussian fluctuations super* 
imposed. The general equation for the ramp pdf is: 

i-o-r-rf] 

Cb7) 

where the incomplete Gamma function is defined as 


i ( 68 ) 

There are 4 constants in equation (22) which are adjust- 
able: a' ,^,c, and cr* La* is an implicit normalizing factor.) 

Antonia et.al. have shown that if a'»1.8 , 8«1.25, c«l,0, and 
0*»O.25, good correlation with experimental results is obtained. 
(For a fuller discussion of the ramp pdf and its properties, see 
reference (1).) .The equation to evaluate the expected values 
(which leads to the moments) for the ramp pdf is : 



(z"**' + (N-iV) 


The effect of changes in the four adjustable constants on 
the skewness and flatness factor is shown in Figures 31 ~ 34. 

It can be seen that the skewness and flatness factor are very 
sensitive to changes in a*. Varying a* from 0.1 to 0.5 will 
yield values of F between 1.0 and 3.7. This sensitivity requires 
that’CT* be selected in accordance with the best-available 

/ p \ 

experimental data. Antonia ' indicates that a value of a*s:0.25 
yields very good agreement with his experimental data. 


with cf* heid constant atg*«0.25 and the other constants 
varied I th^ flatness factor ranges approximately between 2,0 

and 2.5, This is in agreement with Fiedler who predicted 

« 

that a value of 2 was characteristic of the sawtooth-like 
ramp pdf. The range of the skewness was -1,57 to -1.25. 

•The skewness and flatness factor distributions for the 
beta pdf are shows in Figures 27 - 30. The flatness factor 
ranges from 1*66 to 8.64, Since the beta pdf is continuous, 
every interjnediate value can be obtained by the proper selection 
of t and This means that the beta pdf will likewise generate 

the values 2.0-2, 5 for the flatness factor. However, the beta 
pdf generates values of skewness between -2.66 and +2,66. When 

the flatness factor lies between 2*0 and 2.5, the skewness varies 

« 

between -0.82 and +0.82. Recalling that the ramp pdf has skew-* 
ness values between -1.57 and -1.25 for the same range of flatness 
factors, implies that the skewness can be used as the secondary 
criteri^a for distinguishing between the two pdf's. 

The procedure is as follows; if the flatness factor lies 
outside the range 2.0 - 2.5, then always choose the beta pdf 
since the ramp pdf does not produce these values. When the 
flatness is between 2.0 and 2.5, both pdf's can take on these 
values; hence, the skewness must be consulted. If the skewness 
is greater than -0.82, choose the beta pdf; if less than -1.25, 

t 

select the ramp pdf. 

A simpler, though less general, secondary criteron is to 
determine if the skewness is positive or negative. If it is 
positive, choose the beta pdf since the ramp pdf can only produce 
na^ative values for the skewness. If it is negative, the method 
indicated in the last paragraph must be used. 


. The skewness alone cannot be used as a criterion because 
the values generated by the two j/df s are not distinct, but 
overlap. The full range of the skewness values for the ramp 
pdf is -1.57 to -1.25; and the range of the Beta pdf is -2.66 
to +2.66. 

In order to determine values for the skewness and flatness 
factor in a turbulent, reacting flow, a transport equation is 
being developed that will enable the determination of ^quantities. 
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List of Symbols 

pre-exponential constant 
temperature exponent 
instantaneous concentration 
initial concentration 
expected (mean) value of x 
denotes a fuel species 

joint profitability density function of (x,y) 

function of the two-dimensional continuous 
random variable (x^y) 

mean value of g(x#yl 

function of the continuous random variable, x 
Arrhenius reaction rate constant 
constants dependent upon Tmin and Tmax 

laminar Arrhenius reaction rate constant 

mean turbulent Arrhenius reaction rate 
constant 

constant, defined in equation C30) 

denotes maximum value 

denotes minimum value 

order of reaction constant 

denotes oxidizer species 

denotes a product species 

probability density function of x 

joint probability density function of (x,y) 

dimensionless concentration and constant defined 
by equation C30) 

mean dimensionless concentration 

mean square fluctuation of dimensionless 
concentration 

dimensionless temperature 

mean dimensionless temperature 

mean square fluctuation of dimensionless 
temperature 

temperature 
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activation temp®i:atut® 

V6«l. vatianoQ of ix 

V’ value of intogtal defined by equation (44) 

; reaction rate 

I 

laminar reaction rate 

Wi, mean turbulent reaction rate 

i# 


2 reaction rate constant amplification ratio 

2' reaction rate amplification ratio 

f gamma function 

mean of a continuous random variable , x 

VK Kth moment of a random variable 

2 

cf variance of a continuous random variable 

p correlation coefficient 
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AS»PENDIX: Moments oj 
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f pit) Ji ( 

o 
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where 




Integrate Eg. (1) by parts 


Ut pC-t) , 'f-'t 

CX^ *•■ XJ:)‘ pCt)i-i ^ dv^d't 
' jpCi)A-t, l.i. 


(A.fX.f 


I 

- ( y,-t-jya)<ii - fz x^i‘‘.p(i) 


I 

^K*y>* I'v) _ y I .., 3.x, ft"- /;. «)cl-£ =. / 


Use Eq. (3) to eliminate the 

I 




term. 


/ 2 ^ + 4 *■ ) . (■£ ) «/•! *• i£ ' 

O • 


f t t 



I 



16) 


Substitute Eq . C6) into Eq . C5 ) to obtain : 

V.t - X ^ » I 


C7) 


Equation C7) is one equation that contains only A_ X, Xi. 

O JL f /» 

plus constants. 

1+ X, t + 'i (•{; 


X. 


A left J 

* 4 4 




^ X,4 Kw 


C7a) 


To obtain another equation, start with Eq. (2) ; 

I 




a, 4+ 


c/t - *4 
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since , we would like the derivative of the 

exponent to precede the exponential te; 3 . This is done 
by multiplying and adding appropriate terms until 


is obtained, 

f 




di » I 


eft + \ X/C. • dt S 

“ . 0 

o 

I 

J + x,') e ' * dt c ^ X^i -I- X. 

o •'I 

- X 4- X J 
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» 

X^4 X| 4‘ X^ 1 . . 

e - ^ X.-fc + X. 


- X, X X^t + X 
c « 
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(8*) 


C8) 


Egs. C7) and (8) are two equations in three unknowns. The 
third equation may be any one of the Egs. (1), C2) , or (3). 
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Figure 4. Variation of the reaction rate constant 
amplification ratio with the temperature 
fluctuation 
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Figure 5. Variation of the reaction rate constant 
amplification ratio with the temperature 
fluctuation 
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Figure 7. Variation of the reaction rate constant 
amplification ratio with the temperature 
fluctuation 















Figure 12. Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the most-likely bivariate pdf and the 
modified bivariate Gaussian pdf 
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Figure 13. 


t 

Vux'iation of the reaction rate amplification 
ratio with temperature and species fluctuations 
for the most-likely bivariate pdf and the 
modified bivariate Gaussian pdf 






Figure 14 . Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the most-likely bivariate pdf and the 
modified bivariate Gaussian pdf 
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Pigure 15. Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the most-likely bivariate pdf and the 
modified bivariate Gaussian pdf 
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Figure 16. 


Variation of the reaction rate araplification 
ratio with temperature and species fluctuation 
for the most“likely bivariate pdf and the 
modified bivariate Gaussian pdf 
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Figure 17. Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the most-likely bivariate pdf and the 
modified bivariate Gaussian pdf 
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Figure 18. Variation of the reaction rate amplification 
ratio with the maximum temperature for the 
most-likely bivariate pdf and the modified 
bivariate Gaussian pdf 
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Figure 19 . 


Variation of the reaction rate amplification 
ratio with the correlation coefficient for the 
most-likely bivariate pdf and the modified 
bivariate Gaussian pdf 
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Figure 21. Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the beta pdf, m.ost-likely bivariate pdf, 
and the temperature-derived most-likely pdf 
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Figure 22. 


Variation of the reaction rate amplification 
ratio with temperature and species fluctuations 
for the beta pdf, most-likely bivariate pdf, 
and the temperature-derived most-likely pdf 
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V<iriation of the reaction rate amplification 
ratio vdth temperature and species fluctuations 
for the beta pdf, most-likely bivariate pdf, 
and the temperature-derived most-likely pdf 
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Figure 24. 


Variation of the reaction rate amplification 
ratio with temperature and species fluctuations 
for the beta pdf, most-likely bivariate pdf, 
and the temperature-derived most-likely pdf 
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Figure 26. Variation of the reaction rate amplification 

ratio with temperature and species fluctuations 
for the beta pdf, the most-likely bivariate pdf, 
and the temperature-derived most-likely pdf 
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